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ABSTRACT 

Q\ I Just as the ordinary Doppler effect serves as a tool to measure radial velocities of 

. celestial objects, so can the relativistic Doppler effect be implemented to measure a 

CO ■ combination of radial and transverse velocities by using recent improvements in ob- 

serving techniques. A key element that makes a further use of this combination feasible 
is the periodicity in changes of the orbital velocity direction for the source. Two cases 
are considered: (i) a binary star; and (ii) a solitary star with the planetary compan- 
ion. It is shown that, in case (i), several precision Doppler measurements employing 
• the gas absorption cell technique would determine both the total orbital velocity and 

Oh' the inclination angle of the binary orbit disentangled from the peculiar velocity of 

q , the system. The necessary condition for that is the measured, at least with a modest 

■ precision, proper motion and distance to the system. 

CO 

C3 
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1. INTRODUCTION 2 2. THE METHOD 



Measuring transverse (tangential) velocities has always been 
one of the most important and, at the same time, the least 
advanced issues in astronomy. Available methods include 
measurement of either annual parallaxes or proper motions, 
which is possible only for nearby objects. Furthermore, the 
method of proper motions requires knowledge of the distance 
to the object. This Letter shows that, just as the ordinary 
Doppler effect serves as a tool to measure radial velocities of 
celestial objects, the relativistic Doppler effect can be imple- 
mented to measure a combination of radial and transverse 
velocities by using the rapidly improving technique of pre- 
cision Doppler measurements. This combination, as I show 
here, makes it possible to derive the inclination angle of a 
binary star provided that both the distance to the binary 
and its proper motion are measurable. 



The frequency of a spectral line measured at the Earth and 
corrected for the Earth's rotation is supposed to be associ- 
ated with the geocentric coordinate system. It is necessary 
to introduce at least two more reference frames: the frame of 
the source barycentre, S', and the frame of the Solar system 
barycentre, S" . Let the motion of the source relative to 5" 
be described by a (generally variable) velocity = v/c, S' 
moves relative to S" with a (constant in both direction and 
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time) velocity b = V/c, and the Earth moves relative to S" 2.1 2.1. Case (i): a binary 



with a velocity /?©. The relativistic Doppler equation that 



relates the observed frequency of a spectral line, v, to the 



emitted frequency, uo, is given by successive application of 



the usual 4-wavevector transformation (Landau & Lifshitz 



1951) between all the reference frames involved: 



v _ v v 1 v" _ y/l-P 2 Vl~¥ 1 - /3, 



■® cos ( 



U V 1 V" IS 1 — /3r 1 — b r ^/T 



BC£, (1) 



where B absorbs the source's internal velocities (3 and (3 r , 



C absorbs the motion of the source barycentre with velocity 



b — V /c, £ absorbs the Earth's velocity and 9 is the 



angle in S" between /?® and the direction to the source. 



isolation of different types of velocity strictly within the cor- 



responding multipliers, i.e. factorization - makes it impos- 



sible to determine the total and radial velocities separately. 



However, special circumstances, such as periodic changes of 



the direction of the source's velocity, improve the situation. 



Below, I consider two particular cases: (i) the source as a 



binary star and (ii) the source as a star with the planetary 



companion. 



Let us consider a binary in a circular orbit with inclination 
angle i. In this case, the change of the radial velocity (3 r is 
periodic while the total velocity (3 =\ (3 | is constant. Then 
B changes in a periodic fashion as well. Let f3 r reach its 
maximum, /3 r ,max, at an instant t\ so that B has a maxi- 
mum £> max at that instant too. In half a period (the instant 
£2), (3 r = — /3 r ,max so that B reaches its minimum 

^min • 

Ev- 
idently, 



B^ 1 

'-'max 1 '-'mm 



(2) 



i.e. the dependence on /?r,max is cancelled out and the result 
depends on f3 2 only. Once [3 2 is found, /3 r ,max can be derived 
from 



B~ L -B~ l 

*-^mm '-'max 



2(3 r 



It is easy to see that 



At first sight, the very structure of equation (1) - an /3r,max /3sini. 



Therefore equations (2) and (4) yield 

^max ^min 1 



Bmax ~\~ Brain (3 



where 



(3 



( 2B m ax*-' mi ii \ 
mm / 



(3) 



(4) 



(5) 



(6) 



Bmax + B 

Thus both the total velocity (3 and the inclination angle i 
are derivable from the precision Doppler measurements . 

In practical implementation of the above approach, we 
assume that the measured frequency i> is corrected for an- 
nual motion of Earth (term £). Motion of the visible star 
in the binary (with variable velocity (3) and motion of the 



binary barycentre (with constant velocity 6) result in a (vari- 
able) total Doppler shift 

The instants when vq/v reaches its maximum and minimum 
enable us to determine the instant when f3 r = 0. At that 
instant, observations bring us an approximate value of b r 
(by neglecting the higher order terms \b 2 + \p 2 ). Measure- 
ments of the binary's proper motion and distance yield the 
transverse velocity of the binary bt and thus enable one to 
evaluate b 2 — b 2 + b 2 . The difference between vq/v (mea- 
sured at the instant when (3 r = 0) and (l — b r + |& 2 ) yields 
an approximate value of i/3 2 . Obviously, /3 2 is kept constant 
in orbit. By applying the iteration method, the procedure is 
repeated until it converges to give accurate values of b r . As 
soon as b r , b 2 , and /3 2 are found, the measured change of 
vo/v in orbit yields \(3 r \, the modulus of f3 r - Finally, knowl- 
edge of both \/3 r \ — /3sini and (3 2 allows us to evaluate sin i. 
However, if (3 <C b, which would correspond to wide binaries 
or a small mass of the companion, practical implementation 
of the above approach is very difficult, as illustrated by case 
(ii). 
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2.2 2.2. Case (ii): a single star with a planetary 
companion 

Since the motion of a star that has a planetary companion 
orbiting around it is only slightly disturbed by the planet, 
this case is reduced to case (i) of a binary with /3 <C b. In the 
frame S" associated with the barycenter of the Solar Sys- 
tem b is the speed of the source, b r is its radial component, 
and b t is the transverse component. Although f3 <C b, the 
disentangling of /3 r and b r does not make any problem as 
long as the terms linear in f3 and b are only involved. How- 
ever, in this case, as opposite to case (i), determinarion of 
the inclination angle of the planetary orbit turns out to be 
in practice very difficult. The reason is that inaccuracy in 
measuring bt, the transverse component of b, is currently too 
large. To see this, it is worth discussing the issue of accuracy 
in more detail. 

To account for the contribution of the second-order 
terms that enter the expansion of the relativistic Doppler 
formula into a series in v/c, precision of Doppler measure- 
ments of spectral line shifts (translated into the precision of 
radial velocity measurements) should be no less than 

r 1 V 

ov = —v— 

2 c 

= 67 m/s if v = 200 km/s (halo stars), 

= 17 m/s if v = 100 km/s (bulge stars), 

= 1.5 m/s if v = 30 km/s (disk stars). (8) 
Remarkably, accounting for the relativistic Doppler effect 

appears to be within the limits of modern astronomical tech- 
niques (see Sec. 3), even when the measured velocity does 
not exceed 30 km/s. 
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In contrast to the radial velocity, the tangential veloc- 



ity currently cannot be measured with a similar precision. 



However, since vt enters the relativistic Doppler equation 



only as a second and higher order terms, accuracy in mea- 



suring proper motion and distance to the star are allowed 



to be rather modest. For instance, if an inaccuracy in the 



inferred v t does not exceed, say, 10%, as a product of com- 



bined uncertainties in measuring the distance to and proper 



motion of the star, a similar inaccuracy is translated into 



the value of total velocity, provided that the radial veloc- 



ity measuremenets are as accurate as to obey equation (8). 



Although this would still make possible to evaluate the incli- 



nation angle of a binary as described by case (i) , it fails for a 



planetary companion because the current inaccuracy in mea- 



suring b t exceeds the value of (3 that has to be derived. For 



instance, even if SVt = 1 km s 1 , the planetary companion 



induces v typically not exceeding 0.1 km s 1 . Nonetheless, 



since the "signal" (the value of v t ) changes periodically, and 



the "noise" (the value of Vt) is kept constant, an extraction 



of that signal from the noise seems to be although a difficult 



but in principle solvable problem. 



3 3. PRECISION VELOCITY 
MEASUREMENTS 

As can be seen from equation (8), the proposed method re- 
quires high precision in velocity measurements as well as 
stability of the instrumental system. Any seasonal instru- 
mental effects must be carefully eliminated. An idea of the 
appropriate reference frame for measuring velocities, namely 
passing the starlight through an absorbing gas prior to its 
entrance into the spectrograph, was proposed long ago (Grif- 
fin & Griffin 1973) but only recently has been developed 
to a level sufficient to provide high enough precision. The 
modern technique employs a stabilized gaseous iodine (I2) 
absorption cell as a wavelength standard (Libbrecht 1988; 
Marcy & Butler 1992; Cochran & Hatzes 1994; Kiirster et 
al. 1994). The I 2 cell used in the wavelength region 5000-6000 

o 

A, in combination with a high resolving power 1Z provides 
a velocity precision proportional to 7?. _3 / 2 (S/N) _1 , where 
S/N is the signal-to-noise ratio (Hatzes & Cochran 1994). 
The long-term velocity precision achievable with the cur- 
rent techniques is 4-7 m/s (Kiirster et al. 1994) to 3 m/s 
or even better (Butler et al. 1996). For discussion of con- 
ditions that must be controlled in order to achieve such a 
high precision (changes in dispersion and the spectrograph 
point spread function, careful determination of the topocen- 
tric velocity relative to the Solar system barycentre, etc.), 
see e.g. Marcy & Butler (1992). Other effects resulting in 
radial velocity changes, which have been revealed and stud- 
ied during programmes of searches for extra-solar planets, 
include the variability of K giants (e. g. Hatzes & Cochran 
1994), variability among non-Cepheid stars in the instabil- 
ity strip (Butler 1992), and rapid oscillations of Ap stars 
(Hatzes & Kiirster 1994). 

Recently, the stabilized gas absorption cell technique 
has been implemented by several groups (Marcy & But- 
ler 1992; Cochran & Hatzes 1994; Kiirster et al. 1994) in a 
search for Jupiter-like planets. This procedure is commonly 
named 'precision radial velocity measurements', although, 
as is clear from equation (1), it is the changes in the total 
velocity (i. e. including the transverse velocity) that are ac- 
tually being measured as soon as the accuracy of Doppler 
measurements reaches the level given by equation (8). 



4 4. DISCUSSION 

A central point of this letter is that if the accuracy of pre- 
cision Doppler measurements is as high as to comply with 
equation (8), a combination (7) of linear and quadratic in 
velocity terms, and not just radial velocity alone, is actually 
measured. If the radial velocity changes periodically (the 
source is a binary star), the binary parameters, including 
the inclination angle of the orbit, can in principle be dis- 
entangled. In practice, this task might be not easy. For in- 
stance, in an extreme case of a contact binary consisting of 
solar-mass stars (a ~ 1.5 R©, P ~ 0.2 d ), the velocity vec- 
tor changes in At sec seconds by Av = (v 2 /a) At ~ 0.14At scc 
m/s, i.e. by 1.4 m/s every 10 sec, which requires a very short 
exposure time to make measurements. 

Even if the source is a solitary star, precision Doppler 
measurements separated in time can bring important dy- 
namical information, especially if the star is located in a 
dense stellar field such as a globular cluster. Indeed, differ- 
entiation of equation (1) gives the frequency of a spectral 
line (corrected for annual motion of Earth) measured with 
time interval At to be shifted by 



Av / Vn V t 2 VV \ . 
v \ c ca tr / 



(9) 



where V is acceleration and d is distance to the source. A 
similar equation is used in pulsar timing (e.g. Blandford et 
al. 1993), where it is written in terms of pulsar period and 
its derivative. In a globular cluster with the core radius r c — 
0.1 pc and star density within the core p c — 3 • 10 6 M & /pc 3 , 
the dynamical acceleration in the core is V = 4/37rGp c r c = 
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1.7 • 10 _s cm s~ 2 . By neglecting the 2nd and 3rd terms in 
equation (9), one finds that in At yl years the star's velocity 
would change by AV = 5.4Ai yr m/s, i.e. by 5 m/s in a year 
of precision Doppler measurements. Therefore, they seem 
appropriate as a tool to measure gravitational field in the 
cores of globular clusters - an opportunity that so far seemed 
only possible for a few millisecond pusars. 

The prospects for application of the I2 absorption cell 
technique (Sec. 3) seem to be very promising. Combining 
the echelle spectrograph with the I2 cell allows for routine 
observations of stars of magnitude V = 13 (Marcy & Butler 
1992). If the effects of all systematic errors can be reduced 
so that the velocity precision is limited only to the photon 
statistics, then the velocity precision improves as N^ 1 ^ 2 , 
N being the number of detected photons. Fourier analysis 
makes it possible to extract the embedded periodicities with 
amplitudes of about 1 m/s (Brown 1990). This, as is seen 
from equation (8), enables one to measure the relativistic 
correction to the ordinary Doppler effect even when peculiar 
velocities of stars are ~ 30 km/s or less (disk stars). 

Summarizing, this Letter demonstrates in principle a 
possibility to disentangle, for binary stars, the first and sec- 
ond order velocity terms of the relativistic Doppler effect 
and, as a result, to extract the inclination angle of the or- 
bit (by neglecting gravitational fields). This could be done 
by making use of recent substantial improvements in the 
technique of precision Doppler measurements combined with 
proper motion data. However, much further work is needed 
to relax the underlying explicit and implicit assumptions 
and to evaluate the associated corrections as well. 
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